Influence of static Jahn- Teller distortion on the magnetic excitation spectrum of 
PrC^: A synchrotron x-ray and neutron inelastic scattering study 
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A synchrotron x-ray diffraction study of the crystallographic structure of PrCh in the Jahn- 
Teller distorted phase is reported. The distortion of the oxygen sublattice, which was previously 
ambiguous, is shown to be a chiral structure in which neighbouring oxygen chains have opposite 
chiralities. A temperature dependent study of the magnetic excitation spectrum, probed by neutron 
inelastic scattering, is also reported. Changes in the energies and relative intensities of the crystal 
field transitions provide an insight into the interplay between the static and dynamic Jahn- Teller 
effects. 

PACS numbers: 61.10.Nz, 71.7D.Ch, 71.70.Ej, 78.70.Nx 



I. INTRODUCTION 

In recent years there has been strong interest in the 
interplay between electric and magnetic degrees of free- 
dom in lanthanide and actinide dioxides, e.g. UO2, NpC>2 
and Pr02-— These compounds share many similar struc- 
tural and magnetic characteristics, but also present a 
wealth of unusual individual behaviour. For example, 
the longstanding mystery over the nature of the single 
phase transition at 25.5 K in NpC<2 brought about sugges- 
tions first of octupolar ordering, 2,3 then of triakontadipo- 
lar ordering^ Also, a recent study of UO2 by resonant x- 
ray magnetic scattering has revealed electric quadrupolar 
ordering that coincides with a Jahn- Teller distortion and 
the onset of antiferromagnetic ordering. Pr02 exhibits 
an unusually large static Jahn- Teller distortion at 120 
Kj£ as well as a strong dynamic Jahn- Teller effect due 
to the high degree of orbital degeneracy in the ground 
state. 7 In all three compounds magnetoelastic interac- 
tions play a role in determining the low temperature 
properties. However, due to the complexity of such inter- 
actions, the extent of their contribution is not yet fully 
understood^&i^iii 2 . 

In this paper we focus on PrC>2, which exhibits strong 
correlations between its structure and magnetism. At 
room temperature it is paramagnetic and possesses the 
fluorite crystal structure. Below T D — 120 K it under- 
goes a cooperative Jahn- Teller distortion which doubles 
the unit cell along one axis. 6 Magnetic ordering does not 
accompany this transition, but below Xb it is possible to 
influence the population of structural domains through 
the application of a magnetic fields Below T N = 13-5 
K the Pr spins order antiferromagnetically, adopting a 
magnetic structure with two components: one with the 



same unit cell as the fluorite crystallographic structure, 
and the other with the same unit cell as the distorted 
structure below Td Below Tn the populations of both 
structural and magnetic domains can be influenced by 
the application of a magnetic fieldfi 3 - 

We present a synchrotron x-ray diffraction study of the 
crystallographic structure, which distinguishes between 
two previously ambiguous possibilities for the distortion 
below Td . Both possibilities consist of an internal distor- 
tion of the oxygen sublattice, but one involves a shear and 
the other a chiral deformation. We also present tempera- 
ture dependent measurements of the magnetic excitation 
spectrum, probed by neutron inelastic scattering, which 
show that the dynamic Jahn- Teller effect persists over a 
wide range of temperature from 1.8 K to room tempera- 
ture. Changes to the crystal field excitations observed on 
cooling through the Jahn- Teller transition are interpreted 
using a point-charge model of the crystalline electric field 
at the Pr site, drawing on the results of the synchrotron 
x-ray diffraction study. The implications of our results 
for existing models of the dynamic Jahn- Teller effect^ 
are discussed. 



II. X-RAY DIFFRACTION EXPERIMENT 

The fluorite structure (space group Fm3m) exhibited 
by PrC>2 at room temperature is illustrated in Fig. [lja). 
In a recent neutron diffraction study we found that a co- 
operative Jahn- Teller transition at Td = 120 K caused an 
internal distortion of the oxygen sublattice which doubled 
the unit cell along the [100] axis. 6 However, in that study 
the exact configuration of the oxygen ions remained am- 
biguous. Two possible structures were suggested: an or- 
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thorhombic structure (space group Imcb), shown in Fig. 
[IJb), in which the cubic oxygen sublattice was sheared 
along a single crystallographic direction, and a tetrago- 
nal structure [space group 74(1) / acd\, shown in Fig.QJc), 
in which the oxygen ions were displaced chirally, with 
neighbouring chains possessing opposite chiralities. 




(b) 




FIG. 1: (Color online), (a) Pr02 possesses the fluorite crystal 
structure above To- The spheres are the oxygen ions. The 
Pr ions are not shown here, but they occupy positions at the 
corners and face centres of the dashed cube, (b) Sheared 
structure, (c) Chiral structure, (b) and (c) are alternative 
models for the structural distortion that occurs below To- 
The arrows show the directions in which the oxygen ions are 
displaced. 



In the paper published simultaneously with this one, 
Jensen describes a theoretical model which shows that 
the chiral structure is the more stable of the twoJ^ Al- 
though both structures exhibit identical structure factors 
for the Bragg reflections measured previously by neutron 
diffraction, a separate set of reflections exists, which dis- 
tinguishes the two. This set is present for the chiral 
structure and absent for the sheared structure. How- 
ever, the structure factors are extremely weak compared 
to those of the parent fluorite structure. As large sin- 
gle crystals are not available, we decided to exploit the 
high dynamic range achievable with synchrotron x-ray 
diffraction to study these weak reflections. 

We present the results of an experiment performed at 
the European Synchrotron Radiation Facility, Grenoble, 
France on the instrument ID20. The sample used for the 
experiment was a single crystal of mass ~ 1 mg (the same 
sample as used for previous neutron diffraction studies), 
possessing a single flat faceJ^ The sample was enclosed in 
a closed-cycle refrigerator, with a base temperature of 12 
K, which was mounted within a five-circle diffractometer 
with a vertical scattering plane. An incident energy of 
10 keV was used, corresponding to a wavelength of 1.24 

A. 

We studied three sets of Bragg reflections, two of which 
were previously studied by neutron diffraction. The first 



set, corresponding to the fluorite structure, obeys the se- 
lection rule h, k, I all even or all odd. We label this set 
F. These reflections contain scattering intensity from 
both the Pr and O atoms, but the cross section for x- 
rays is dominated by the contribution from Pr due to its 
much larger atomic number. Since the Jahn- Teller (JT) 
distortion is predominantly an internal distortion of the 
oxygen atoms, the fluorite reflections remain strong be- 
low Td- The second set, corresponding to the distorted 
structure below Td , and arising purely from oxygen scat- 
tering, obeys the selection rule h = {In + l)/2, k = odd, 
I = even, I ^ (or permutations thereof). We label 
this set JT(I). This set is common to both the sheared 
and chiral structures, so it cannot help us distinguish be- 
tween the two. The third set, which also arises purely 
from oxygen scattering, was not previously studied by 
neutron diffraction. This set obeys the selection rule h 
= odd, k = 2n, I = An (or permutations thereof), and is 
present for the chiral structure only. We label it JT(II). 

First, we selected three Bragg reflections typical of 
each of the sets described above, and measured the tem- 
perature dependence of these using #-scansi£ over a range 
of temperatures between 20 K and 150 K. A Ge (111) 
analyzer crystal was used for these scans in order to op- 
timize the resolution and dynamic range. The integrated 
intensity and the full width at half maximum (FWHM), 
obtained by fitting a Lorentzian profile to each #-scan, 
are plotted as a function of temperature for the three 
peaks in Fig. [5] Above Td only the fluorite peak (422) 
is present. Below Td, the oxygen peaks (43|) and (630) 
appear. The presence of the JT(II) peak (630) immedi- 
ately rules out the sheared structure and provides strong 
evidence in favour of the chiral structure. 

Next, we measured all the accessible peaks in each set 
using 0-scans at 20 K to determine the crystallographic 
structure factors. These scans were done without the 
Ge (111) analyzer crystal in order to improve the count- 
ing statistics. The integrated intensities were obtained 
by fitting Lorentzian profiles and were subsequently cor- 
rected for the geometric Lorentz factor and the absorp- 
tion factor. The peaks in set F were also corrected for 
attenuation occurring from an aluminium plate of 1 mm 
thickness that was placed in the beam to prevent detector 
saturation. Finally, the intensities were multiplied by a 
universal scaling factor to allow direct a comparison with 
the squared structure factors |T(Q)| 2 calculated for the 
chiral structure. The scaling factor was chosen such that 
the intensity of the (630) reflection matched its squared 
structure factor exactly. 

Table U compares the corrected scaled intensities of the 
observed peaks with the calculated squared structure fac- 
tors |T(Q)| 2 for the chiral structure, assuming each oxy- 
gen ion is displaced by 0.0726 A in the direction shown 
in Fig.njc)j& The observed intensities agree with the cal- 
culated |T(Q)| 2 to within a factor of 4 for sets JT(I) and 
JT(II), and to within an order of magnitude for set F. 
Better agreement cannot be expected due to the difficulty 
in estimating the absorption factor away from the specu- 
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FIG. 2: (Color online). Temperature dependence ol the in- 
tegrated intensity (blue circles) and FWHM (red squares) ol 
typical peaks from sets F, JT(I) and JT(II). The insets 
show t9-scans below and above To (open circles: T — 20 K, 
closed circles T > 120 K). Integrated intensities are normal- 
ized to the incident beam intensity. The insets show detector 
counts, which are not normalized. 



lar condition. There is also an additional uncertainty on 
the measured intensity of the peaks in set F due to the 
possible nonuniformity of the thickness of the aluminium 
attenuator. In addition to the measurements shown in 
Table |l] we checked several reflections whose indices did 
not satisfy the selection rules for sets F and JT(II). All 
of these were absent, as expected. 



III. NEUTRON SCATTERING EXPERIMENT 

In a previous experiment we probed the low tempera- 
ture excitation spectrum of PrC>2 using neutron inelastic 
scattering on a powder sampled The spectrum showed 
two main features: sharp peaks characteristic of crystal 
field transitions at energies above 100 meV, and a broad 
band of scattering in the range 10-100 meV which we 
interpreted as magnetic scattering from vibronic excita- 



Set 


Reflection 


|F(Q)| 2 (fm 2 ) 


Intensity (fm 2 ) 




(311) 


2.6 x 10 4 


2.0 x 10 4 




(222) 


1.7 x 10 4 


1.5 x 10 5 


F 


(333) 


1.7 x 10 4 


1.6 x 10 5 


(422) 


2.4 x 10 4 


1.3 x 10 5 




(402) 


1.4 x 10 4 


5.3 x 10 4 




(531) 


1.4 x 10 4 


4.1 x 10 4 




(2|1) 


4.9 


2.5 






5.5 


2.8 


JT(I) 


(63|) 


4.3 


13.6 




(43§) 


3.3 


5.9 




(6|1) 


5.3 


15.0 




(630) 


0.3 


0.3 


JT(II) 


(401) 
(201) 


0.2 
0.04 


0.2 
0.03 




(432) 


0.05 


0.19 



TABLE I: Comparison between the squared structure factors 
F(Q)| 2 predicted for the chiral structure [Fig. [ljc)] and the 
integrated intensities of the observed x-ray reflections at T = 
20 K. F denotes fluorite reflections, JT(I) denotes reflections 
from the oxygen displacements that are common to both the 
sheared and chiral structures, and JT(II) denotes reflections 
arising solely from the chiral structure. Uncertainties are not 
given for the integrated intensities, as these are dominated by 
the uncertainty in the absorption factor, which is difficult to 
quantify. 



tions associated with the dynamic Jahn- Teller effect. We 
developed a simple model of the magnetoelastic coupling 
between a single phonon mode and a crystal field exci- 
tation, which yielded a qualitative agreement with our 
results. However, at the time we were unaware of the 
existence of a static Jahn- Teller distortion, and therefore 
chose a phonon mode consistent with cubic symmetry 
for our model. Other researchers have also used cubic 
symmetry as the basis for more complex multiphonon 
models^ In the light of our present knowledge, that the 
tetragonal distortion involves only an internal displace- 
ment of the oxygen atoms, leaving the cubic Pr lattice 
undisturbed, the choice of cubic symmetry is not unrea- 
sonable. However, the distortion lowers the crystal field 
symmetry at the Pr site and the large magnitude of the 
oxygen displacements is expected to cause a significant 
splitting of the cubic Tg ground state, as well as an alter- 
ing of the positions and intensities of the higher energy 
crystal field transitions. An examination of the temper- 
ature dependence of the excitation spectrum is therefore 
in order. 

We report the results of a neutron inelastic scatter- 
ing experiment performed at the ISIS spallation source, 
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Didcot, UK on the high energy transfer (HET) time-of- 
flight chopper spectrometer. A powder sample of PrC>2 
was prepared from a starting material of PrgOn, which 
was baked at 1000 °C for several hours, then annealed 
in flowing oxygen at 280 °C for ~ 5 days. The prod- 
ucts were checked by x-ray diffraction and no trace of 
residual PrgOn could be found. A sample of mass 25.83 
g was sealed in an aluminium foil packet and mounted 
inside a temperature-controlled top-loading closed-cycle 
refrigerator. The spectrum was measured at 7 K, 50 K, 
80 K, 100 K, 110 K, 115 K, 120 K, 130 K, 165 K and 
200 K, using two different incident energies at each tem- 
perature: Ei = 80 meV (chopper frequency 350 Hz) and 
Ej = 250 meV (chopper frequency 500 Hz). An iden- 
tical set of measurements (although at a reduced set of 
temperatures: 7 K, 80 K, 120 K and 165 K) was made 
on a sample of CeC>2 of mass 25.98 g, sealed in a simi- 
lar foil packet, to estimate the non-magnetic background 
due to phonons and multiple scattering in PrC>2 ^ Vana- 
dium spectra were used to calibrate the detectors and 
convert the scattering intensity into absolute units of 
cross-section, i.e. mbsr~ 1 meV _1 (Pr ion) -1 . 

Figure [3Ja) shows PrC>2 spectra measured in the low 
angle detector banks with an incident energy of Ei = 250 
meV, corrected for absorption and self-shielding. The 
single peak is due to a crystal field transition from the 
ground state to the excited state. Above Tb the 
peak is centred at 121 meV, but as the sample is cooled 
through the Jahn- Teller transition the width of the peak 
decreases sharply and its centre moves rapidly to higher 
energies, reaching 132 meV at 7 K. The data are fitted us- 
ing a Lorentzian profile on a background estimated from 
the CeC>2 data. The fit parameters are used to plot the 
temperature dependence of the peak centre, FWHM and 
integrated intensity in Fig.[3](b). 

The change in energy of the peak below Td can be un- 
derstood qualitatively by constructing a simple model of 
the crystalline electric field, ignoring the spatial extent 
of the electron wavefunctions and treating the Pr 4+ and 
2 ~ ions as point charges. This type of model underes- 
timates the separation of the crystal field levels because 
it disregards the overlap between the electron orbitals. 
However, it gives a reasonable idea of the relative energy 
separation of the levels and also their relative spectral 
weights. If the separation of a pair of levels is known 
experimentally, the point-charge model can be scaled to 
match, thus enabling us to predict the approximate po- 
sitions of other levels. 

Figure HJa) shows a point-charge simulation of the ex- 
citation spectrum for the fluorite phase of PrC>2 above 
Td- The simulation takes into account the resolution 
function of the HET spectrometer at an incident energy 
of 250 meV and assumes a sample temperature of T = 
K, but the intensity is not normalized to absolute units 
of cross section. In the fluorite phase, the ground state 
is a Tg quartet and the excited state is a Ti doublet. 
The separation of these has been scaled to match the ob- 
served transition energy of 121 meV. Figure \Mjo) shows 
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FIG. 3: (Color online), (a) Crystal field transition from the 
ground state to the IV excited state in Pr02, measured by 
neutron spectroscopy. The red circles are data from the low 
angle detector banks of the HET spectrometer. Data taken at 
different temperatures have been offset vertically for clarity. 
The solid black line is a fit to the data using a Lorentzian pro- 
file on a background estimated from the Ce02 data (dashed 
blue line), (b), (c) and (d) show the temperature dependence 
of the peak centre, the FWHM and the integrated intensity, 
respectively. 



the simulated spectrum from a point charge model of the 
chiral phase below Td, using the same scaling factor as 
for the fluorite phase. It is immediately obvious that the 
T7 level has moved to a higher energy, in agreement with 
our experimental observations (although the magnitude 
of the observed shift on cooling through Td is rather less 
than the point-charge prediction). The model also shows 
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FIG. 4: (Color online), (a) Point-charge model for the fluo- 
rite crystal structure above Td. The model has been scaled to 
match the observed Ts —* transition energy of 121 meV. 
The thick black vertical lines represent the peak intensities. 
The thin blue line is a simulation of the neutron scattering 
excitation spectrum, but is not in absolute units of cross sec- 
tion, (b) Point charge model for the chiral structure below 
Td, using the same scaling factor as the fluorite model. 



that the spectral weight of the FV level decreases in the 
chiral phase, which accounts for the discrepancy between 
model and experiment in our previous study^ Lastly, the 
Tg ground state is predicted to split into two doublets, 
separated by ~ 21 meV. 

The results of the point-charge model lead us to ex- 
amine the low energy spectrum of PrC>2, to see if we can 
detect the splitting of the T 8 ground state below Td . Fig- 
ure [5] shows the difference between the PrC>2 and CeC>2 
spectra measured in the low angle detector banks with an 
incident energy of Ei = 80 meV. CeC>2 data for temper- 
atures intermediate between actual measurements were 
estimated by interpolation. The broad band of scatter- 
ing above 10 meV is clearly present both below and above 
Td ■ We observe that the scattering above 35 meV is rel- 
atively independent of temperature. However, the region 
below 35 meV is strongly temperature dependent with 
a maximum that shifts from ~ 28 meV at 7 K to ~ 20 
meV at 100 K, before becoming quasielastic close to Ti> 
The energy of the maximum at 7 K is in good agree- 
ment with the prediction of the point-charge model for 
the chiral phase and also with other prediction s 14 ^ 5 for 
the splitting of the r§ ground state in the distorted phase 
below Td- This provides strong evidence that the broad 
band of scattering contains a component arising from a 
crystal field transition between the split doublets of the 
r 8 ground state. 

In order to extract the temperature dependence of 
the splitting in an unbiased way, we constructed a line- 
shape to fit the broad band of scattering from three sym- 
metrized Gaussians and a symmetrized Lorentzian, as 
shown in Fig. [TJ] Peak A represents quasielastic scatter- 
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FIG. 5: (Color online), (a)-(d) Pr0 2 spectra with the Ce0 2 
background subtracted (circles) , showing the variation of the 
broad band of vibronic scattering with temperature. The 
black line is a fit, constructed from four separate peaks (see 
Fig. |6j , which is used to determine the temperature depen- 
dence of the crystal field transition between the split doublets 
of the Tg ground state in the chiral phase. The centre of the 
transition, determined from the fit, is indicated in each plot 
by a blue arrow. 



ing from the ground state, peak B represents the tran- 
sition between the two doublets of the split T% ground 
state, and peak C represents a broad continuum of vi- 
bronic scattering. Peak L (the symmetrized Lorentzian) 
was added to improve the quality of the fit at low en- 
ergies. Each peak was weighted by a detailed balance 
factor to take into account the thermal population of the 
levelsJ£ In order to determine the width and amplitude 
of peak B, the lineshape was initially fitted to the spec- 
trum at 7 K. We then fitted the lineshape simultaneously 
to all the spectra for temperatures between 7 K and 130 
K, keeping the width and amplitude of peak B fixed. The 
other peak widths, centres and amplitudes were free pa- 
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rameters, but were not allowed to vary with temperature. 
The only parameters allowed to vary with temperature 
were the centre of peak B and the amplitude of peak L. 




50 100 150 

Temperature (K) 

FIG. 6: (Color online). Temperature dependence of the 
ground state splitting in the chiral phase. The circles are 
the results of the fit to all spectra between 7 K and f30 K, 
where the centre of peak B was allowed to vary with temper- 
ature. The dashed line is a guide to the eye. The inset shows 
the functions used to construct the lineshape which was fitted 
to the data. A, B and C are symmetrized Gaussians, whereas 
L is a symmetrized Lorentzian. All are weighted by a detailed 
balance factor. 



confirms that the displacement of the oxygen atoms in 
the Jahn- Teller distorted phase is chiral and a neutron 
inelastic scattering experiment which probes the crystal 
field transitions above and below the Jahn- Teller transi- 
tion temperature Tb . The cubic T% ground state is found 
to split into two doublets in the distorted phase, and the 
splitting is found to be ~ 20-30 meV, in good agreement 
with a simple point charge model of the crystal field. 

Our results indicate that existing models of the dy- 
namic Jahn- Teller effect^ should be revised to include 
the effect of the static distortion. The accompanying pa- 
per by Jensen^ presents such a model, which exhibits 
an improved agreement between the calculated and mea- 
sured intensities of the crystal field transitions. Our neu- 
tron scattering data show that vibronic excitations due to 
the dynamic Jahn- Teller effect are largely independent of 
temperature, despite the presence of a static Jahn- Teller 
transition, ft is likely that this is because the magnetoe- 
lastic energyiii is comparable to the crystal field split- 
ting of the ground state in the distorted phase. Hence, 
the static Jahn- Teller effect does not quench the dynamic 
Jahn- Teller fluctuations. 



The solid black lines in Fig. [5] show the fitted line- 
shape, which agrees well with the data. Figure [5] shows 
the temperature dependence of peak B. The peak clearly 
moves to lower energies as the temperature rises, becom- 
ing quasielastic near Tq. 

IV. SUMMARY AND CONCLUSIONS 

We have presented the results of two experiments on 
Pr02: a synchrotron x-ray diffraction experiment which 
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The flat face of the crystal was close to the (111) plane, but 
did not coincide directly with a plane of high symmetry. 
Therefore, a Laue image was used to align the crystal with 
respect to the instrument. 

A (9-scan in a five-circle vertical diffractometer is the equiv- 
alent of an oj-scan in a four-circle horizontal diffractometer. 



This type of scan involves the rotation of the crystal about 
an axis normal to the scattering plane, while keeping the 
scattering angle 26 fixed. 

CeC>2 possesses the fluorite structure and has a similar 
lattice parameter to Pr02- The neutron scattering lengths 
of Ce and Pr differ by no more than 10%, but the Ce 4+ ion 
has no 4/-electrons, making it non-magnetic. Hence, Ce02 
provides a good estimate of the non-magnetic scattering 
from Pr02- 



